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Abstract. Carbon burning is a crucial process for a number of important astrophysical scenarios. The lowest 
measured energy is around E cm =2.1 MeV, only partially overlapping with the energy range of astrophysical 
interest. The currently adopted reaction rates are based on an extrapolation which is highly uncertain because 
of potential resonances existing in the unmeasured energy range and the complication of the effective nuclear 
potential. By comparing the cross sections of the three carbon isotope fusion reactions, 12 C+ 12 C, 12 C+' 3 C and 
13 C+ 13 C, we have established an upper limit on the molecular resonance strengths in 12 C+ 12 C fusion reaction. 
The preliminary results are presented and the impact on nuclear astrophysics is discussed. 



1 Introduction 

In 1960 Almqvist, Kuehner and Bromley discovered sev- 
eral resonances in collisions between 12 C nuclei that pro- 
vided the first evidence for the existence of nuclear molecules. 
For at least three energies, E c m =5.68, 6.00 and 6.32 MeV, 
they observed resonances in the yield of collisional by- 
products: protons, gamma radiation, alpha particles and 
neutrons. The resonances had characteristic widths of about 
120 keV, which corresponds to a lifetime of 5xl0~ 21 sec- 
onds, a value more than 10 times greater than the typ- 
ical time span for a simple collision ! 1121 . In the follow- 
ing years, as the measurements were pushed down towards 
lower energies, the discoveries of such resonances contin- 
ued as far as the lowest measured energy. For instance, the 
most recent published measurement of the 12 C+ 12 C fusion 
reported a strong resonance at E cm =2.14 MeV. However, 
the reported width (less than 12 keV) is much less than the 
common 100 keV observed at higher energies|3|. 
Carbon burning driven by the 12 C+ 12 C fusion, is a crucial 
process for a number of important astrophysical scenarios, 
such as the formation of white dwarfs, nucleosynthesis in 
massive stars, carbon ignition in type la supernovae and in 
superbursts[4|. The important energy range spans from 1 
MeV to 3 MeV in the center of mass frame, which is only 
partially covered by the measurements of the 12 C+ 12 C fu- 
sion reaction. The extrapolation is the only resource avail- 
able to obtain the fusion reaction rate for the astrophysi- 
cal applications. The current adopted 12 C+ 12 C fusion reac- 
tion rate is established based on the modified S factor(S*), 
which is defined as, 

S*(E) = o-(E)Ee M i +046£ . (1) 



An averaged S* factor of 3xl0 15 MeV*b was obtained by 
fitting the data measured by Patterson[5], Spinka|6| and 
Becker[7|. This averaged value was extrapolated down to 
the energies of astrophysical interest with a simple assump- 
tion that the averaged modified S factor is constant at the 
sub-barrier energies (8). 

The strong and relatively narrow resonances continue down 
to the lowest energies for which measurements have been 
made. To present, there is nothing known about the exact 
energies and strengths of these resonances in the energy 
region which is significant in the astrophysical context. 
Therefore, the currently adopted carbon fusion rate based 
on the averaged S* factor is highly uncertain. 
Superbursts are long, energetic, rare explosions in Low 
Mass X-Ray Bursts. These bursts are considered to be trig- 
gered by the unstable carbon burning in the ash left over 
from the rp-process on the surface of a neutron star. The 
ash is heated up with the heat generated in the crust of a 
neutron star and ignites the carbon in the ash. However, 
with the currently adopted carbon fusion rate, superbust 
models fail to explain the ignition of carbon at the column 
depth inferred from observation[9]. Another difficulty is 
that the amount of 12 C in the ash does not seem to be suf- 
ficient to trigger the superburst. Inspired by the strong res- 
onance observed at E c m .=2.14 MeV, a hypothetical reso- 
nance with a similar strength was proposed to enhance the 
fusion reaction rate at Tg- 0.5 K with a factor of over two 
orders of magnitude and thereby alleviate the discrepancy 
between superburst models and observations IfTOll . 
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Fig. 1. The experimental S* factors of three carbon isotope fusion 
reactions, 12 C+ 12 C I3I5I7I1 II . 12 C+ l3 C fl2lT5l and 13 C+' 3 CfT3l. 
The averaged S* factor of 12 C+ 12 C(S*=3xl0 16 MeVb), shown 
as red line, is recommended based on the data of Becker, Patter- 
son and Spinka (not shown here). There are some data with large 
error bars in Spillane's data. In this paper, only those with rela- 
tive errors less than 70% are shown. The hypothetical resonance 
proposed by Copper et al. is also shown. 



2 Empirical relationship among the carbon 
isotope fusion reactions 

In contrast to the striking resonances in the 12 C+ 12 C fu- 
sion reaction, other carbon isotope fusion reactions, such 
as 12 C+ I3 C and 13 C+ 13 C, behave more regularly. Little res- 



onance feature has been observed in these two systems 1 121131 . 
The shape of the fusion cross sections at sub-barrier ener- 
gies is primarily dominated by the Coulomb barrier pene- 
tration effect. To remove this effect and reveal more details 
of nuclear interaction, the cross sections of all three carbon 
fusion systems are converted into the S* factors using eq[T] 
The advantage of using this conversion over the traditional 
astrophysical S factor is that the cross section ratios among 
the three systems are preserved in this approach. 
The experimental S* factor for all three carbon isotope fu- 
sion systems are shown in Fig[T] There are several impor- 
tant features in these carbon isotope fusion systems: 

1. The 12 C+ 12 C cross sections are bounded from above 
by the cross sections of the other two carbon isotope 
fusion systems. 

2. All the maxima of the 12 C+ 12 C fusion cross section 
(E,.=3.08, 4.28, 4.92, 5.67, 5.96 and 6.26 MeV) match 
remarkably well with the relatively smooth cross sec- 
tions observed for the other two isotope combinations, 
especially 13 C+ I3 C. 

3. Overall, the 12 C+ 13 C cross sections are the highest among 
the three carbon isotope fusion systems. The deviation 
between 12 C+ I3 C and I3 C+ 13 C occurs in the energy of 
3.5 to 5 MeV (excluding the 4.25 MeV resonance in 
12 C+ 13 C) and is less than 30%. 

4. The 4.25 MeV resonance in 12 C+ 12 C is also present in 
12 C+ 13 C. 

The Nogami-Imanishi model has suggested the cou- 
pling effect to the 12 C(2+,4.44 MeV) is crucial to the for- 



mation of 12 C+ 12 C molecular resonances! 14]. In the 12 C+ 13 C 
and 13 C+ 13 C fusion reactions, the valence neutron of 13 C 
introduces an additional freedom for excitation besides the 
12 C core excitation. Despite this complication, the maxi- 
mum deviation among the 12 C+ 12 C peak cross sections, 
12 C+ 13 C and 13 C are less than 30%. Based on this strong 
correlation, we claim: 

1. All major resonances have the same reaction mecha- 
nism; 

2. This reaction mechanism should be very similar to the 
process which fuses 12 C+ 13 C and 13 C+ I3 C. 

3. The effect arising from the valence neutron in 13 C is 
not strong compared to this reaction mechanism. 

The Imanishi calculation may provide a qualitative expla- 
nation to the 12 C+ 12 C experimental data. However, no quan- 
titative description has been achieved yet. 



3 The 12 C+ 13 C experiment at Notre Dame 
and test of extrapolating models for heavy 
ion fusion reactions at extreme sub-barrier 
energies 

Up to now, measurements of 12 C+ 12 C have reached down 
to 2.1 MeV at which the cross section is below 1 nb. How- 
ever, measurements of 12 C+ 13 C and 13 C+ 13 C stop at en- 
ergies above 3 MeV. The lowest measured cross sections 
are at the level of jjb. It would be interesting to extend the 
measurements of 12 C+ 13 C and 13 C+ 13 C towards lower en- 
ergies and examine the empirical correlations among the 
three systems. 

The correlations among the three carbon isotope systems 
provide a possibility to establish an upper limit for the im- 
portant 12 C+ 12 C. If all the major resonances in 12 C+ 12 C, 
including those to be discovered at lower energies, are driven 
by the same reaction mechanism, one can expect to estab- 
lish an extrapolation based on the resonant cross sections 
observed at higher energies. To achieve a more reliable ex- 
trapolation within the energy range of 1 to 3 MeV, it is 
necessary to push the measurements of the carbon isotope 
fusion reactions towards lower energies. 
As the first step, we have studied the 12 C+ 13 C fusion within 
the center of mass energy range of 2.6 MeV to 4.8 MeV 
using the FN tandem accelerator at Notre Dame. In this 
experiment, the I2 C( 13 C, 24 Na)p channel was studied via 
decay spectroscopy. The total fusion cross section was ob- 
tained by correcting the measured cross section using the 
branching ratio from Ref. lfT2l . The lowest cross section 
reached in this measurement is about 20 nb, a factor of 
50 smaller than which Dayras achieved lfT2l . The experi- 
mental detail can be found in Ref.[15|. The new I2 C+ 13 C 
data are shown together with the Dayras measurement in 
fig. [T] and [2] In general, the shape of the new data agrees 
with the Dayras measurement except in the range of 3.3 
to 3.5 MeV where the new data is a little lower than the 
Dayras measurement. It may be explained as a fluctuation 
in the branching ratio incurred by the nuclear structure in 
the compound nucleus, 25 Mg. It is also observed in fig. [I] 
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Fig. 2. The experimental S* factors for 12 C+ 13 C I12I15I . The 
ESW fit to the Dayras data (blue points) is shown as blue line 
while the Akyiiz-Winther calculation is shown as magenta line. 
The ESW prediction match the new measurement studied at 
Notre Dame (green points). 
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Fig. 3. The experimental S* factors for 12 C+ 12 C I3I5I7I11I . The 
ESW model is fit to the major peak cross sections at E r =3.08, 
4.28, 4.92, 5.67, 5.96 and 6.26 MeV. The ESW fit is shown as blue 
line. The Akyiiz-Winther calculation is also shown as magenta 
line. It is clear that the resonant cross section at 1.5 and 2.14 
MeV are too high to be realistic. 



that the new 12 C+ I3 C data still lie above the 12 C+ 12 C cross 
sections measured by Spillane et al. within the overlapped 
energy range. 

The fusion systems with smooth cross sections, includ- 
ing I2 C+ 13 C and 13 C+ 13 C, can be reasonably described by 
coupled-channel calculations which include the couplings 
between the relative motion of the colliding nuclei and sev- 
eral nuclear intrinsic motions as well as transfer reactions. 
One challenge in this approach is the determination of the 
effective nuclear potential. So far, there is no universal po- 
tential which may provide reliable extrapolation down to 
the extreme sub-barrier energies. For example, the cou- 
pled channel calculation with the Akyiiz-Winther poten- 
tial, a global potential with a Woods-Saxon shape, provides 
a good description of the fusion cross sections around the 
Coulomb barrier. However, it was found in the hindrance 
study that, at lower energies, the calculation always over- 
predicts the fusion cross sections lfT6l . Right now, there is a 
ongoing effort trying to model all three carbon isotope fu- 
sion reactions using a unified double folding potential ifTTl . 
In contrast to the standard approach, the Equivalent Square 
Well (ESW) model has been chosen traditionally to model 
the sub-barrier fusion cross sections and provide extrap- 
olation within the energy ranges of astrophysical interest. 
Although the physical justification of the model may not 
be as clear as the coupled channel calculations used in 
the sub-barrier fusion study, this model has achieved fits 
to measured data that are superior to those obtained from 
most other models l[T8l . To examine the predictive power of 
the ESW model for the carbon isotope fusion reactions, the 
ESW extrapolation determined in an earlier measurement 
at higher energies lfT2ll is shown in Fig. [2] together with the 
new 12 C+ 13 C data. The excellent agreement proves the su- 
perior predictive power of the ESW model. As a compari- 
son, a couple channel calculation using the Akyiiz-Winther 
potential was carried out with the code, CCFULL[ 19|. The 
13 C(3/2~, 3.684 MeV) with strong B(E2) is included in the 
calculation. The result is shown in Fig. [2] The calculation 



significantly over predicts the fusion cross section at the 
energies less than 4 MeV. For example, the fusion cross 
section at 2.6 MeV is over predicted by a factor of 3.4. 

4 Upper limit on the 12 C+ 12 C molecular 
resonance strengths 

The resonance strength is directly related to the peak cross 
section since all the nuclear molecular states are expected 
to have a similar width. The correlation between the 12 C+ 12 C 
peak cross sections and the cross sections of the other two 
carbon isotope fusion systems provides a possibility to model 
the resonance strengths. The ESW model has been applied 
to fit those major peak cross sections ( E r =3.08, 4.28, 4.92, 
5.67, 5.96 and 6.26 MeV) which match the 12 C+ 13 C and 
13 C+ 13 C fusion cross sections. The ESW extrapolation is 
shown in Fig. [3] together with the experimental data. Even 
though the fit is based on the peak cross sections above 3 
MeV, the extrapolation to lower energies provide an upper 
bound on most of Spillane's data well, except the strong 
resonance observed at 2.14 MeV. The resonant cross sec- 
tion at this energy is about a factor of 40 higher than the 
ESW extrapolation. The preliminary result from the 12 C+ 12 C 
experiment [20 1 at Naples supports our finding. 
Another, even safer, upper limit can be estimated by coupled- 
channels calculation that are based on the he Akyiiz-Winther 
potential. This approach is known to over-predict the fu- 
sion cross sections at extreme sub-barrier energies. The 
calculation is carried out using CCULL[19| by including 
the coupling to the 12 C(2 + ,4.44 MeV). The result is also 
shown in Fig. [3] The two resonances at 1.5 MeV and 2.14 
MeV are about a factor of 10 higher than the Akyiiz-Winther 
calculation. 

We cannot rule out the possibility of a strong reso- 
nance in the unmeasured energy range, formed by a mech- 
anism different than the ones we have discussed here. The 
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established upper limit cannot account for such a reso- 
nance. However, there is no such evidence thus far. In prin- 
ciple, there are some sharp resonances existing in the 24 Mg 
which decay by gamma emission. It is possible that those 
resonant cross sections are much higher than the ESW ex- 
trapolation, but the resonance strengths are too small to 
affect the 12 C+ 12 C fusion reaction. 



5 Impact to nuclear astrophysics 

Even though we do not know the exact location of the po- 
tential resonance in the unexplored energy range, with the 
established upper limit, we can claim the real 12 C+ I2 C fu- 
sion rate is much less than Cooper et al. proposed. Taking 
the over-predicted S* factors calculated with the Akyiiz- 
Winther potential, the real rate should not be more than a 
factor of 6 stronger than the currently adopted fusion rate 
based on the S* factor in the energy range of 1 to 3 MeV. 
The upper limit established with the ESW fit, the best fit 
we have achieved, is only 20% higher than the currently 
adopted rate. The averaged S* factor of Spillane's data 
tends to suggest a smaller average than the current value. 
Therefore, the more than a factor of 40 enhancement as 
proposed in Ref.[ 10] is impossible. 

The upper limit posts an additional constraint on the cur- 
rent superburst models. Since the carbon burning rate can 
not be as high as Cooper et al. proposed, there must be 
some unknown physics processes not being included in the 
superburst models. For example, an additional heating pro- 
cess in the neutron star crust could raise the peak tempera- 
ture in the ash to meet the required ignition condition. It is 
also possible that some process other than carbon burning 
triggers the superburst. 

6 Summary and outlook 

Carbon fusion is an important process. The molecular reso- 
nance complicates the prediction of the carbon fusion rate 
at the energies of astrophysical relevance. By comparing 
the cross sections of three carbon isotope fusion reactions, 
12 C+ 12 C, 12 C+ 13 C and 13 C+ 13 C, an empirical relationship 
among the three systems is found. With this relationship 
and the Equivalent Square Well (ESW) model, for the first 
time, the resonant cross sections in the 12 C+ 12 C fusion re- 
actions is quantitatively described. Consequently, an upper 
limit on the molecular resonance strengths in 12 C+ 12 C fu- 
sion reaction is established. This upper limit serves an ad- 
ditional constraint for the superburst models. 
To provide a more accurate prediction at extreme sub-barrier 
energies, it is urgent to improve the experimental technique 
so that we can push the measurements of the carbon iso- 
tope fusion reactions towards lower energies. Right now, 
the Notre Dame group is building a 5 MV single end ac- 
celerator with an ECR source to provide more than 40 p//A 
carbon beam. We are also improving detection techniques 
to increase efficiency as well as selectivity [21 1. These ap- 
proaches are expected to provide better experimental data 
in the near future. The proposed underground laboratory 



with a high current heavy ion accelerator will be the ideal 
place to provide the ultimate answer to the carbon fusion 
rate. 
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